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Abstract 


A  numerical  scheme  for  the  approximation  of  a  parameter-dependent  problem 
is  said  to  exhibit  locking  if  the  accuracy  of  the  approximations  deteriorates 
as  the  parameter  tends  to  a  limiting  value.  A  robust  numerical  scheme  for  the 
problem  is  one  that  is  essentially  uniformly  convergent  for  all  values  of  the 
parameter.  We  develop  precise  mathematical  definitions  for  these  terms,  give 
their  quantitative  characterization  and  prove  some  general  theorems  involving 
locking  and  robustness.  A  model  problem  involving  heat  transfer  is  analyzed 
in  detail  using  this  mathematical  framework  and  various  related  computational 
results  are  described.  Applications  of  our  theory  to  some  different  problems 
involving  locking  are  presented. 
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1 .  Introduction. 


The  mathematical  formulation  of  various  problems  involves  dependency  on  a 
crucial  parameter  t  arising  out  of  physical  considerations.  For  example, 
plate  and  shell  models  involve  the  thickness  d  of  the  plate  or  shell,  the 
analysis  of  elastic  materials  in  general  involves  the  Poisson  ratio  v,  heat 
transfer  involves  the  ratio  of  conductivities  p  in  different  directions, 
etc. 

The  numerical  approximation  of  such  parameter-dependent  problems  may  suf¬ 
fer  when  t  lies  close  to  a  limiting  value  tQ  (thickness  d— >0,  Poisson 
ratio  v — >0.5,  ratio  p  of  conductivities  is  very  small  (p — >0),  etc.). 

Most  a  priori  error  estimates  yield  optimal  asymptotic  convergence  rates  for 
t  >  tQ  fixed.  Since  these  estimates  are  not  uniform  in  t,  a  degeneration 
often  occurs  for  values  of  t  close  to  t^.  This  is  manifested  in  actual 
computations,  where  the  error  may  not  decrease  at  the  predicted  rate  for  t 
close  to  tQ  for  most  practical  choices  of  the  discretization  parameter.  We 
refer  to  this  phenomenon  as  locking.  Various  examples  of  locking  have  been 
reported  in  the  engineering  and  mathematical  literature,  see  for  example  [1], 
[6]  -  [10],  [14]  and  others. 

Locking  involves  the  "shifting"  of  the  asymptotic  range  of  the  calcu¬ 
lations  and  will  eventually  disappear  when  the  level  of  discretization  is 
increased  enough,  depending  upon  the  strength  of  the  degeneracy.  Unfortunate¬ 
ly,  using  brute  force  may  lead  to  an  infeasible  level  of  discretization 
required  before  convergence  can  be  observed.  Moreover,  the  prediction  of  the 
discretization  level  for  a  required  error  tolerance  is  now  quite  complicated. 

The  most  ideal  remedy  to  locking  is  to  employ  a  method  that  is  robust . 
l.e.,  one  which  is  more  or  less  uniformly  convergent  for  all  t.  Various 
robust  methods  (particularly  mixed  methods)  have  been  proposed  and  analyzed  in 


the  context  of  locking,  see  for  example  [7]  for  the  plate  problem,  [1],  [9] 
for  a  beam  problem,  [6],  [14]  for  Poisson  ratio  locking,  among  others. 


Generally,  in  the  above  mentioned  papers,  locking  was  addressed  in  an 
ad  hoc  manner.  In  this  paper  our  primary  goal  is  to  develop  a  systematic 
mathematical  approach  which  will  allow  precise  characterization  of  locking  and 
robustness  of  a  method  as  well  as  a  quantitative  measure  of  the  locking  and 
robustness  strength.  In  the  next  section  we  present  the  definitions  of  these 
concepts. 

Our  definitions  are  quite  general  in  nature,  so  that  one  can  treat 
various  different  types  of  locking  phenomena  from  the  same  point  of  view.  The 
examples  analyzed  in  Sections  2-5  (both  theoretical  and  computational  1 
illustrate  the  need  for  this  flexibility. 

In  the  next  section,  we  also  derive  various  general  theorems  for  a 
special  class  of  locking  problems  which  are  important  in  the  context  of  appli¬ 
cations  and  have  been  studied,  for  example,  in  [6].  These  problems  have  a 
well-defined  associated  limit  problem  for  t— »tg,  satisfying  certain  proper¬ 
ties.  An  example  is  the  case  of  elasticity,  where  we  get  the  Stokes’  problem 
in  the  limit  as  the  Poisson  ratio  v— >0.5.  In  [6],  it  was  shown  that  a 
necessary  condition  for  tUe  absence  of  locking  is  the  satisfaction  of  a  suit¬ 
able  approximabi 1 ity  condition  for  this  limit  problem.  We  show  in  Theorems 
2.2,  2.3  that  this  condition  is  sufficient  as  well,  under  certain  conditions. 

Section  3  contains  detailed  results  on  locking  and  robustness  for  the 
model  problem  Introduced  in  Section  2.  Section  4  shows  how  the  non-quasi- 
uniformity  of  meshes  used  can  be  interpreted  in  terms  of  locking.  In  Section 
5,  we  present  a  simplified  analysis  for  the  beam  problem  which  was  analyzed  in 
[1],  [9].  In  Section  6,  we  discuss  a  similar  simplified  analysis  that  is  pos¬ 
sible  for  the  Stokes’  problem,  detailed  results  for  which  (using  the  theorems 
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from  Section  2)  are  presented  in  [5]. 


2.  A  model  problem  and  the  definitions. 

To  show  the  main  features  of  locking,  we  first  discuss  an  illustrative 
model  elliptic  problem  which  arises  in  the  study  of  (for  example)  heat  trans¬ 
fer  through  highly  orthotropic  materials  with  the  conductivity  being  very 
different  in  two  perpendicular  directions.  Let  the  conductivities  in  the  x^ 
and  x2  directions  be  k^  =  1  and  k^  =  1/t,  respectively.  We  will  be 
particularly  interested  in  the  case  that  t  is  close  to  0. 

We  assume  homogeneous  Dirichlet  boundary  data  on  and  Neumann  data  on 

fN  (with  r  =  rD  u  rN).  The  corresponding  partial  differential  equation  we 


study  is 

(2.  1) 

2  2 
a  u  i  a  u 

f(x,y)  in 

Q 

(2.2a) 

u  =  0 

on  rD 

(2.2b) 

au 

an  =  g 

c 

°n  rN. 

where  n 

c 

is  the  usual 

conormal . 

Let 

yn) 

=  (u  e 

H^fl)  |u  =  0  on 

Td>  .  Then 

this  problem  may  be  put 

the  variational 

form: 

Find  u  €  H^j(fi) 

such  that 

for  all  v  e  H^(fl) , 

(2.3) 

bJ(u, v) 

=  F(v) , 

where 


F(v)  = 


fvdx^dx^  + 

n 


gvds. 


Let  n  =  (-1,1)  for  simplicity.  We  will  be  particularly  interested  in 
the  case  that  =  0,  f  =  0,  and  g  is  chosen  so  that  the  true  solution 


(unique  up  to  a  constant)  is  given  by 
(2.5) 


1  ,  “VTx2 

u^  =  sin  x^e 


As  t— >0,  it  is  seen  that  uj  has  a  well  defined  limit  sin  x^  which  is 
constant  in  the  xg  direction. 

We  also  consider  (2.3)  with  B*  replaced  by  its  rotated  form 


(2.6)  B^(u,v) 


fau  +au  ] 

fav  av  ] 

ifau 

au  ] 

fav 

-3v  11 

laxi  +  ax2J 

K  + 

t[ax2 

“axj 

K 

a*Jj 

which  corresponds  to  the  case  where  the  directions  of  the  orthotropy  do  not 
coincide  with  the  x^,x2  By  choosing  g  appropriately,  we  obtain  as 

the  solution  of  our  new  variational  problem  the  rotated  version  of  (2.5) 


(2.7) 


ufc  =  sin(x1+x2)e 


-v't(x2— Xj) 


2 

which  is  once  again  unique  up  to  a  constant.  As  t— >0,  ufc  becomes  constant 
in  the  t j  =  x2  -  xj  direction. 

We  now  approximate  the  above  problems  by  the  finite  element  method. 

First  we  consider  the  h-version  on  a  uniform  square  mesh  on  fl  with  1,  4 

and  16  square  elements.  In  Figures  2.1  -  2.3  we  have  plotted  the  percentage 

relative  error  versus  number  of  degrees  of  freedom  in  a  log-log  scale,  for 

1  2 

p  =  1,  2  and  4.  In  each  case,  results  are  shown  for  both  u^.  and  ut>  for 

_  1  —  Q  1 

t  =  10  and  10  .  The  dotted  lines  show  the  error  in  the  H  norm,  while 


the  solid  lines  represent  the  error  in  the  energy  norm  -  their  behavior  is 


HKMTIVE  E11UOK 


very  similar  (see  Theorem  2.4). 
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Figure  2. 1.  Relative  error  in  and  for  the  h-version, 
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1  2 

Figure  2.2.  Relative  error  in  u  and  u.  for  the  h-version 


Figure  2.3.  Relative  error  in  and  u  for  the  h-version,  p  =  4. 

It  is  observed  that  for  t=  10  ,  both  cases  lead  to  the  error  decreas- 

_0 

ing  at  the  predicted  asymptotic  optimal  rate.  However,  when  t  =  10  ,  only 

the  unrotated  case  shows  the  same  behavior.  The  error  for  the  rotated  case 
hardly  decreases  -  this  is  described  as  a  locking  effect. 


Figure  2.4. 


and  u  '  for  the  p-version  using  one  square  element 


Relative  error  in  uj 


- NUMBER  OF  DEGREES  OF  FREEDOM  — *• 

Figure  2.5. 

1  2 

Relative  error  in  ut  and  for  the  p-version  using  four  elements. 

In  Figures  2.4  -  2.5  the  p-version  has  been  used,  with  h  =  2,1  (i.e., 

1  and  4  square  elements).  Figure  2.4  shows  that  the  p-version  with  one 

1  2 

square  is  free  of  locking  for  both  u^  and  u^.  since  the  value  of  t  has  no 

appreciable  effect  on  the  observed  rate  of  convergence.  In  Figure  2.5,  it  is 

2  “6 

seen  that  the  error  curve  for  u^  with  t  =  10  is  shifted  upwards  compared 

to  the  other  curves,  although  the  slope  of  this  curve  remains  the  same.  This 

suggests  that  the  order  of  convergence  remains  the  same,  but  the  constant  in 

the  error  estimate  is  larger  for  t  small.  (This  "shift"  will,  of  course, 

occur  for  the  h-version  as  well,  but  will  only  be  visible  when  h  is  of  the 

same  order  as  t.  Hence,  in  practice  it  is  not  observed.) 

Figure  2.6  compares  the  h-version  (p  =  1,2)  and  p-version  (one 
2  1 

square)  for  u^.  in  terms  of  the  H  (n)  norm  error.  This  clearly  illustrates 
the  difference  in  terms  of  locking.  It  shows  that  the  p-version  is  more 
robust . 
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- 1  =  10'°  p  =  1,  2  — 

/  h-version 
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- - o  t  =  10  ’  p-version 
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t  =  10  p  =  2 
_  h-version 


t  =  10  p-version 

t=10b  p-version  j  .  i 


—  NUMBER  OF  DEGREES  OF  FREEDOM  — ► 

_  2 

Figure  2.6.  Comparison  of  the  h-  and  p-version  accuracies  for  ut- 

In  Section  3,  we  shall  prove  various  results  to  explain  these  phenomena. 
It  should  be  noted  that  the  above  numerical  results  are  only  valid  for  the 
practical  range  of  h  and  p  we  have  worked  with.  If  refinement  is  con¬ 
tinued,  the  asymptotic  convergence  rate  would  eventually  be  observed  for  any 
fixed  t. 

We  now  develop  mathemat ical ly  precise  definitions  of  locking  and  robust¬ 


ness. 


Assume  first  that  the  parameter  t  is  an  element  of  a  set  S  which  is 
taken  to  be  a  semi-open  interval  of  the  form  (a,b]  (with  the  limiting  value 
of  t  being  a).  For  each  t,  let  B  (•••)  be  a  bounded  symmetric  bilinear 
form  defined  over  a  Hilbert  space  V  such  that  the  energy  norm 

1/2 

(2.8)  IMIe  t  =  (Bt(u,u) ) 

satisfies 


(2.9) 


where 


C1(t )||u||v  <  ||u|lE  t  <  C2(t)||u||v 

0  <  C,,C„  <  oo,  though  C  (t)/C  (t)  could  — or  0  as  t— >a. 

12  \  c* 


8 


(Here,  IMiy  represents  the  usual  norm  of  the  quotient  space  if,  for  example, 
we  have  a  pure  Neumann  problem. ) 

Let  H  =  {H^.}  be  a  family  of  linear  spaces,  c  H  c  V,  where  H  is  a 
compactly  imbedded  subspace  of  V.  For  each  Ht>  let  ||  »||^  denote  an  asso¬ 
ciated  norm  and  let  =  {u  e  H.  ,  ||u||  <  1  > ,  where  ||u||„  S  C||u||.  .  In  our 

t.  t,  t  H  t 

examples,  Ht  may  often  be  identical  to  H. 

Ue  are  interested  in  the  approximation  of  "exact  solutions"  that  lie  in 
Hf  To  this  end,  for  each  t,  let  =  {V^>  be  a  family  of  finite¬ 
dimensional  subspaces  of  V  with  the  dimension  N  being  independent  of  t. 

N 

(We  assume  that  a  V  is  given  for  each  N  e  N,  with  N  being  unbounded.  ) 

N  N 

For  ufc  €  Ht>  we  consider  the  sequence  of  approximations  u^.  e  V  given  by 

(2.10)  Bt(ut.v)  =  Bfc(ut,v)  V  v  e  v£, 

i.e.,  a  sequence  of  approximations  to  the  problem  ¥  associated  with 

bilinear  form  Bt  and  solution  set  Hfc.  The  family  {?t>  can  be  identified 

with  an  extension  procedure  5,  i.e.,  a  rule  which  (for  each  t)  gives  a 

method  regarding  how  the  dimension  N  of  the  approximate  subspace  is 

increased.  It  is  seen  that  (2. 10)  defines  a  projection  of  the  space  into 

N 

Vfc.  However,  the  exact  form  of  (2.10)  may  not  allow  us  to  consider  certain 
problems,  for  example,  those  involving  inhomogeneous  Dirichlet  boundary  condi¬ 
tions.  A  more  general  form  of  (2.10)  is  therefore 

(2.11)  u?  =  BtV 

N  N 

where  Bt  :  H^. — » V  is  a  projection  operator.  The  definitions  that  follow 

are  equally  valid  for  (2.11)  as  well. 

In  case  we  are  using  (2.11),  we  will  assume  the  existence  of  an  energy 
norm  corresponding  to  (2.8)  such  that 
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(2. 12) 


Ut  ~  Ut  E,  t 


<  C  inf 

v€^t 


|Ut-V|l 


E,  t 


for  all  N, t  and  C  Independent  of  t.  Note  that  otherwise  (2.12)  follows 
immediately  from  (2.10). 


Remark .  We  introduced  the  set  by  help  of  the  norm  ||*||^.  This  charac¬ 
terization  is  a  very  special  one.  We  can,  instead,  define  to  be  any 

reasonable  set.  We  will  do  this  in  Theorem  3.5. 

We  denote  =  Sx>V  and  let  R  c  Rq  be  a  subset  of  the  pairs  (t,N) 
under  consideration. 

Our  final  component  is  a  family  §  =  {E^.}  of  measures  for  the  error 
E^.  :  V — »(R.  We  are  interested  in  the  errors  Et(ut~u^)  for  problems 
Let  us  mention  that  some  of  Bt>  H^.,  3^,  Efc  may  be  independent  of  t  (and 
usually  are) . 

For  ary  R,  we  now  define  the  locking  ratio  L( t , N)  =  L(t,N,R, H, ?,g,  rj) 
(0  <  7)  <  b-a)  for  the  problem  {7^} ,  a  <  t  £  b,  by 


(2. 13) 


L(t, N)  = 


IT  f 

su5„  WV 

uteHt  _ 


inf 

(t, N)eR  Ut6Ht 


sup  Et(ut-u^) 

B 


where  R  =  {(t,N),  t  >  a  +  T)>  nR. 


Remark.  The  locking 
value  t  to  the  best 
are  characterized  by 
denominator  over  R 

V 

arbitrarily  close  to 
technicality).  This 
the  problem  is  least 


ratio  compares  the  performance  of  the  method  at  parameter 

possible  performance  for  reasonable  values  of  t,  which 

t  >  a  +  7).  In  most  applications,  the  infimum  in  the 

is  the  same  as  that  over  R,  so  that  tj  can  be  taken 

zero  (and  appears  in  the  definition  essentially  as  a 

N 

is  because  typically,  the  error  in  u^.  is  smallest  when 
singular. 
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Let  us  illustrate  this  ratio  by  considering  the  bilinear  forms  (2.4), 

(2.6)  with  R  =  {10  \  10  ^,10  ^,10  ^}x^V.  Let  (1  *  1,2)  be  spanned  by 

the  single  element  u*,  i  =  1,2  (see  (2.5)  and  (2.7)),  and  let  Et(u)  = 

II UH  E,  t ' 

“■6 

We  take  ?  to  be  the  h  extension  procedure  (and  let  ti  =  10  ).  Then 

1  2 

L  (t,N)  and  L  (t,N)  are  shown  for  various  p  in  Tables  2.1,  2.2  using 
1,  4,  9  and  16  square  elements. 

It  is  observed  that  for  p  =  1,  L*(t,N)  remains  close  to  unity. 

Similar  results  are  observed  for  higher  p  for  L^  (not  reproduced  here). 

2 

On  the  other  hand,  it  is  seen  that  for  all  p,  L  (t,N)  seems  to  be  unbounded 
as  t — >0.  (For  t  fixed,  the  locking  ratio  is  of  course  bounded,  but  can  be 
very  large. ) 

1  2 

Table  2.1.  Locking  ratio  L  (t,N)  and  L  (t,N)  for  the  h-version, 
p  *  1  (17  =  10  )  and  the  energy  norm 
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p  =  2 


p  =  3 


7 

20 

39 

64 

1. 

00 

1. 

00 

1. 

00 

1.00 

1. 

64 

2. 

67 

2. 

92 

2.97 

1. 

77 

4. 

81 

10. 

06 

15.90 

1. 

77 

4. 

87 

10. 

67 

18.76 

1. 

00 

1. 

00 

1. 

00 

1.00 

1. 

04 

1. 

46 

2. 

17 

2.40 

1. 

05 

1. 

79 

4. 

98 

10.43 

1. 

05 

1. 

79 

5. 

11 

11.74 

1. 

15 

1. 

00 

1. 

00 

1.00 

1. 

01 

1. 

73 

1. 

83 

1.50 

1, 

00 

7. 

83 

15. 

17 

15.50 

1. 

00 

9. 

01 

43. 

33 

101.56 

for 

the 

P- 

vers 

ion 

for 

one  and 

elements.  We  see  that  for  one  element,  L  stays  bounded  while  for  16 
elements  it  becomes  unbounded. 
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2  — R 

Table  2.3.  Locking  ratio  L  (t,N)  for  the  p-version  (tj  =  10  ) 

and  the  energy  norm 


#  of 
elem. 

t\N 

3 

7 

11 

16 

22 

10" 1 

1.02 

1.00 

1.00 

1.  15 

1.35 

IQ"2 

1.00 

1.64 

1.04 

1.01 

1.00 

1 

IQ'4 

1.00 

1.77 

1.05 

1.00 

1.00 

0 

1 

CD 

1.00 

1.77 

_ 

1.05 

1.00 

1.00 

#  of 
elem. 

t\N 

24 

64 

104 

160 

10'1 

1.00 

1.00 

1.00 

1.00 

10-2 

1.41 

2.97 

2.39 

1.52 

16 

1 

0 

1.53 

15.91 

10.43 

15.56 

CO 

1 

0 

1.53 

18.76 

11.74 

101.53 

Let  us  now  present  the  definition  of  locking. 


Definition  2. 1.  The  extension  procedure  (F  is  free  from  locking  over  the 
region  R  for  the  family  of  problems  {3*  ,  (t,N)  €  R>  with  respect  to 
solution  sets  K  and  measures  8  iff  there  exists  0  <  77  <  b  -  a  such  that 


(2. 14) 


lim  sup]  sup  L(tf  N,  R,  H,  5,  6,  tj)  1  = 
N-*«o  t  t  J 


=  M  <  00. 


(t, N)eR 


?  shows  locking  of  order  f(N)  (where  lim  f(N)  =  «)  iff  for  some  tj, 

N-»oo 


(2. 15) 


0  <  lim  supj  sup  L(t,N,R,H,  =  C 

N-»eo  ^  t  ' 


<  ». 


(t, N)eR 


For  the  csLse  that  C  is  bounded  (respectively,  infinite),  we  say  that 
the  order  of  locking  is  at  most  (respectively,  at  least )  f(N).  □ 
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Related  to  the  above  definition,  we  may  also  define  robustness  as 
follows. 

Definition  2.2.  The  extension  procedure  ?  is  robust  over  the  region  R  c  RQ 
for  the  family  of  problems  {tP^,  (t,N)  €  R>  with  respect  to  solution 

sets  }f  and  measures  &  iff 

N 

lim  sup  sup  E.  (u  -u.  )  =  0. 

N-*»  t  u  eK?  t  1  l 
(t,N)eR 

It  is  robust  with  uniform  order  g(N)  iff 

(2-16)  UN.«UPif  T  U^H8  Et(Ut'Ut)]irNl)  =  c  <  "• 

(t, N)eR 

where  g(N)— >0  as  N— n 

Definitions  2.1  and  2.2  allow  great  flexibility  in  the  way  various  compo¬ 
nents  may  depend  on  t.  This  is  because  as  illustrated  by  the  examples  in 
this  and  the  succeeding  sections,  different  formulations  require  different  t 
dependencies. 

The  measure  in  our  definition  can  take  various  forms,  for  example, 

the  maximum  stresses,  stress  at  a  point,  etc.  The  two  main  choices  discussed 
here  will  be  the  V  norm  and  the  energy  norm. 

As  t— »a,  the  approximation  may  deteriorate  due  to  reasons  other  than 
locking.  For  example,  the  exact  solution  of  the  Reissner-Mindl in  plate  model 
with  uniform  load  has  a  boundary  layer  of  increasing  strength  as  t — >0.  Hence 
the  deterioration  of  the  approximation  is  due  to  the  loss  in  regularity  of  the 
solution  as  well  as  locking.  The  numerical  resolution  of  these  two  effects  is 
based  on  different  strategies.  Our  definition  permits  us  to  effectively 
isolate  the  locking  effects  by  sufficiently  restricting  the  set  H  to  exclude 
unsmooth  solutions. 
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By  locking,  we  typically  understand  the  behavior  of  the  extension  proce¬ 
dure  when  t  is  much  smaller  than  N  *  (for  a  =  0).  Hence,  instead  of 
considering  the  region  RQ  =  (0,b]xvV  (Figure  2.7a),  we  may  be  interested  in  a 
region  of  the  type  R  =  {(t,N),  0  <  t  ^  bN>  shown  in  Figure  2.7b.  This  is 

why  we  have  R  explicitly  in  our  definition.  For  the  examples  discussed  in 
this  paper,  however,  we  are  able  to  prove  all  our  results  for  the  entire 
region  R  =  RQ. 


X/N — ►  - 1  /  N 


Figure  2.7a.  The  region  R^.  Figure  2.7b.  The  region  R. 

We  will  often  need  to  bound  below  the  denominator  in  (2.13).  In  this 
connection,  the  following  definition  of  the  n-width  will  be  useful.  (Here 
H8  is  the  unit  ball  in  H.  ) 

The  n-width  of  the  set  H  in  V  with  respect  to  the  measure  E  is 
defined  as 

(2.17)  d..( V,  H,  E)  =  infjsup  inf  E(u-v)|M  c  V,  dim  M  =  N> 

N  M  W  veM 


with  M  being  a  linear  manifold.  When  E(u)  is  given  by  ||u||y,  we  denote 
this  by  d..(V,  H) . 

N 

The  choice  of  the  measures  ||»||  and  Efc  are  not  independent  if  we  want 
to  study  problems  which  maJce  sense  mathematically.  We  will  restrict  our 
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attention  to  choices  of  H,g,R  such  that  for  any  (t,N)  €  R, 


(2. 18)  dN(V,Ht,Et)  <  CFa(N). 

where  FA(N)— >0  as  N— »a>,  F^  independent  of  t,  C  Independent  of  N 

and  t. 

Condition  (2.18)  ensures  that  there  are  extension  procedures  for  which  it 
is  at  least  possible  that  the  robustness  can  be  achieved.  Related  to  (2. 18) 
is  the  following  condition,  which  ensures  that  our  choice  of  extension 
procedure  has  a  chance  to  avoid  locking:  For  any  (t,N)  e  R, 

(2.7)  C  F  OOS  sup  inf  Et(ut-V)  SC  F  (N>. 

VHt  veVt 

(t.N)eR 

where  FD(N) — >0  as  N — » oo,  F_  independent  of  t,  C  independent  of  N 

D  D 

and  t.  We  say  that  {Jt,  5.6.R}  is  Fg-admissible  if  (2.19)  is  satisfied. 

Note  that  (2. 19)  can  be  an  unduly  broad  condition  to  assume  -  for  example,  if 
E,  =  II  •  II  —  .  ,  then  (2.19)  already  implies  robustness  with  uniform  order  FD, 
as  well  as  absence  of  locking. 

We  now  formulate 


Assumption  A.  There  exists  0  <  v  £  b  -  a  such  that  for  (t.N)  e  R 

T) 

(2.20)  C  F <N)  S  sup  E  <u  -uj)  <  C  F  (N). 

utsHt 

( t, N)eR 

n 

where  FQ(N)— >0  as  N— >a»,  FQ  independent  of  t,  C^.C^  >  0  Independent 

N 

of  t.N,  but  in  general  depending  on  t>  (ut  is  the  finite  element 
solution).  □ 


Assumption  A  says  that  there  exists  am 
dure  is  robust  with  uniform  order  ^(N). 


R  on  which  the  extension  proce- 
t? 

Under  this  ass’^pf i^n,  the  locking 
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ratio  L  takes  a  simpler  form  and  Definition  2.1  is  equivalent  to  saying  that 
5  is  free  from  locking  iff 

(2.21)  lim  sup  I  sup  sup  E.(u  -u^) (F_(N) )  M  =  M  <  « 

N-»«  ^  t  u  eH®  J 

(t, N)eR  1 

(a  similar  reformulation  holds  for  (2.15)). 

The  following  theorem  relates  the  concepts  of  locking  and  robustness  and 
follows  immediately  from  their  definitions. 

Theorem  2.1.  Let  S,  RQ  =  SxjV,  R  c  Rq,  H,  be  given  and  let  Assumption  A 

hold.  Then  the  extension  procedure  ^  is  free  from  locking  over  the  region 
R  iff  it  is  robust  over  the  region  R  with  uniform  order  Fg(N).  Moreover, 
if  f(N)  is  such  that 

f(N)F  (N)  =  g(N)-»0  as  N — >oo, 

then  5  shows  locking  of  order  f(N)  iff  it  is  robust  with  maximum  uniform 
order  g(N).  □ 

Note  that  5  is  non-robust  iff  it  shows  locking  of  order  (FQ(N))  1.  In 
this  case  we  say  that  ?  shows  complete  locking. 

Let  us  now  discuss  a  special  case  of  the  bilinear  form  B^(*,*)  in  terms 
of  which  many  locking  problems  can  be  formulated.  Let  V  be  imbedded  in  the 
Hilbert  space  W,  with  norm  ||*|lw  and  inner  product  (•,♦).  We  assume,  for 
t  €  S  =  (0, b] 

(2.22)  Bt(u,v)  =  a(u, v)  +  i(Cu,Cv) 

is  a  symmetric  bilinear  form  defined  on  VxV,  where  a( • , • )  is  a  symmetric 
bilinear  form  satisfying 

(2.23)  CjIuHy  S  (a(u,u))1/2  <  C2|u||y 


17 


and  C  .  V — »W  Is  a  linear  form  satisfying 


(2.24)  l|Cv||M  <  C3l|v||v. 

Here  C^.C^.C^  are  constamts  independent  of  t  (a(*,*)  amd  C  could  be 
made  t -dependent  as  long  as  (2.23)  -  (2.24)  were  satisfied).  We  have  chosen 

(a,b]  =  ( 0 , b ]  for  convenience. 

(2.22),  (2.23),  amd  (2.24)  give  immediately  that  for  amy  t  €  (0,b],  the 
following  version  of  (2.9)  holds 

(2.25)  CjHully  <  ||u||Eit  <  C2t“1/2||u||v 

We  shall  give  examples  of  problems  that  satisfy  the  form  (2.22)  in 

sections  3,  5,  6.  For  a  number  of  such  problems  amd  the  set  K  =  {H^.},  Ht  c 

H  c  V,  it  is  useful  to  define  lim  H2  =  H?  c  H  c  V.  We  define  it  as  follows: 

t-0 

B  B 

uQ  €  Hq  iff  there  exists  a  sequence  {ut>,  ut  e  Ht’  e  S,  and  a  constamt 

C  independent  of  t  such  that 

(2. 25)  l|ut!lE,t  -  C 

(2.27)  || u^  -  Uglly— >0  as  t— >0. 

Assume  now  u^  e  H^,  then  by  definition  there  are  ut  e  such  that 

(2.26)  amd  (2.27)  hold.  Hence 

r"Cut»u  s  c 

so  that 

(2.28)  ||  Cut  II  w  =  0(t1/2) 

from  which  it  follows  that  for  amy  u^  €  H^  we  have 

(2.29)  CuQ  =  0. 

In  many  exaimples  it  turns  out  that  Cu^. — >0  at  a  rate  0(t)  instead  of 


(2.28)  (see  the  beam  example  in  Section  5  and  Poisson  locking  in  Section  6), 
i.  e.  , 

(2.30)  l|Cut!lW  =  0(t)  (or  o(t))- 


Obviously  (2.30)  is  satisfied  if  Condition  (a)  below  holds. 


Condition  (a).  Let  ||  *  ||  be  such  that  ||u[|_  < 

H^,  there  is  a  uQ  e  n  (uQ  depending  on 

(2.31)  || ufc  -uQ||t  5  Ct1/2 

with  C  independent  of  t  and  u^. 


C||u||  .  Then  for 
ufc )  such  that 


any 


€ 


□ 


Condition  (a)  is  important  in  the  context  of  various  theorems  and  charac¬ 
terizes  an  aspect  of  the  regularity  of  the  problem. 

We  now  prove  some  general  theorems  for  the  form  (2.22)  which  hold  when 
the  spaces  are  assumed  independent  of  t. 

Theorem  2.2.  (A)  Let  H  =  (Ht>,  Hfc  c  H  c  V,  9  •  {9  } ,  9  =  {vj*},  = 

V**,  g  =  {Et>,  Efc(u)  =  ||u||  y,  S  =  (0,  b] ,  R  =  {  ( t ,  N) ,  0  <  t  <  b^>  c  RQ  and 

let  lim  =  H^.  Assume  further  that  the  quadruple  {X, 9,&, R>  is 
t->00 

FQ-admissible.  Then  the  extension  process  9  is  free  from  locking  over  R 


only  if 


(2.32) 


g(N) 


sup  inf 
ueH®  weVN 
Cw=0 


llu  -  w||v  <  CFQ(N). 


More  generally,  it  shows  locking  of  order  f(N)  only  if 
(2.33)  g(N)  S  CFq ( N ) f ( N ) . 

Moreover,  if 


(2.34) 


g(N)  S  CFq ( N ) f ( N ) 
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then  the  extension  process  shows  locking  of  at  least  order  f(N). 

(B)  Suppose,  in  addition,  that  Condition  (a)  holds.  Then  5  is  free 
from  locking  over  R  iff  (2.32)  holds.  It  shows  locking  of  order  f(N)  iff 

(2.35)  C1FQ(N)f (N)  5  g(N)  S  C2FQ(N)f(N). 


Proof.  First,  it  is  easy  to  see  that  Assumption  A  is  satisfied  with  the  same 
Fq  because  for  t  >  y  the  energy  norm  ||*|lg.  t  is  equivalent  to  the  norm 
||*||y.  Suppose  there  is  no  locking  (i.e.,  ?  is  free  of  locking)  over  R. 
Let  N  be  fixed.  Then  by  (2.21),  for  amy  (t,N)  we  have  for  u^.  e  H“, 


I' ut  ~  ut  1  v  *  CFq(N). 

B  B 

Let  now  uQ  €  HQ  amd  ufc  e  Ht  be  such  that  (2.26)  amd  (2.27)  hold.  Because 


v“ 


N, 


is  finite-dimensional,  there  exists  a  subsequence  of  {ut>  such  that 


N  N  ,,N  , 

u.  — > U-.  €  V  in  V  amd 

t  j  u 


l uo  ”  uo II V  *  CFo(N)* 


with  CuQ  =  CuQ  =  0.  Hence  (2.32)  holds.  (2.33)  follows  similamly  by  replac¬ 
ing  FQ(N)  by  F  (N)f(N). 

Next  suppose  that  the  locking  is  of  order  o(f(N)).  Then  by  (2.33) 


g(N)  <  o(FQ(N)f (N) ) , 


which  contradicts  (2.34).  Hence  the  locking  is  of  at  lea^t  order  f(N). 

Suppose  now  that  Condition  (a)  holds.  We  have  to  show  that  if  (2.32) 
holds  there  is  no  locking.  So  let  for  amy  u^  e  «?•  uo 6  4  n  Hfc  be  such 
that  (2.31)  holds.  Then  we  may  write 

(2.36)  ufc  =  u0+(ut-uQ)  =  u0  +  Zt, 


where  €  H^  satisfies 
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Hence  we  have 


|ut-utlE.t  5  C  l|Ut'V|lE,t 


5  C  |uO**t'u|lE,t 

2  C{  lnfj(  luo"vlllv*  lnf“ 


w/vr"  0 

Cv  =0 


v7" 


t  "  VE.t}’ 


using  the  fact  that  €  HQ  and  hence  Cu^  =  0.  This  gives,  by  (2.32)  and 
(2.25) 


<  C{FQ(N)  +  t*1/2  inf  ||z  -v  ||v>. 
v2€V^ 


Since  satisfies  (2.37),  we  use  the  Fg-admissibi 1 ity  (2.19)  to  get 

-1/2  1/2 

<  C{Fq(N)  ♦  t  X/V Fq(N)> 

<  CFq(N). 

Then  we  see  that 

»ut-\»v  S  C'ut-^«E.t  S  CF0(N> 

so  that  on  R,  L(t,N)  remains  bounded,  i.e.,  there  is  no  locking  in  the  V 
norm. 

Suppose  now  that  (2.35)  holds.  Then  we  may  replace  Fq(N)  in  the  above 
by  FQ ( N ) f ( N )  so  that 

llut  ~  ut(V  5  Cllut-utllE.t  S  CF0lN,f(N) 

which  shows  locking  of  order  at  most  f(N).  Also,  since  (2.34)  is  satisfied. 


the  locking  is  of  order  at  least  f(N),  i.e.,  of  order  f(N). 

Conversely,  if  the  locking  is  of  order  f(N),  then  (2.33)  holds.  If 
g ( N )  =  o(FQ(N)f (N) ) ,  then,  just  like  before,  we  establish  that 

«ut-u”||v  <  o  C  Fq  CN)f(N) ) 

so  that  the  locking  is  of  order  at  most  o(f(N)),  a  contradiction.  □ 

We  have  actually  proved  a  stronger  result,  namely, 

Theorem  2.3.  (A)  Let  the  conditions  of  Theorem  2.2  (A)  hold.  Then  ?  is 

free  from  locking  with  respect  to  the  energy  norm  over  R  only  if  (2.32) 
holds.  It  shows  locking  of  order  f(N)  only  if  (2.33)  holds.  It  shows 
locking  of  order  at  least  f(N)  if  (2.34)  holds. 

(B)  Let  Condition  (a)  be  true.  Then  ?  is  free  from  locking  over  R 
iff  (2.32)  holds.  It  shows  locking  of  order  f(N)  iff  (2.35)  holds. 

Proof.  As  before,  Assumption  A  is  satisfied  so  that  in  both  cases  we  can  use 
the  same  Fq(N).  The  rest  of  the  proof  is  the  same.  □ 

Theorems  2.2  and  2.3  show  that  under  certain  circumstances,  we  need  only 
check  condition  (2.32)  (or  (2.35))  to  determine  whether  or  not  locking  exists. 

The  energy  norm  and  the  V  norm  are  two  important  error  measures.  The 
relation  between  them  in  terms  of  locking  is  expressed  in  the  following 
theorem. 

Theorem  2.4.  (A)  Let  the  assumptions  of  Theorem  2.2  (A)  hold.  Then  5  is 

free  from  locking  over  the  region  R  c  R^  with  respect  to  the  V  norm  if  it 
is  free  from  locking  with  respect  to  the  energy  norm.  It  shows  locking  of 
order  f(N)  in  the  V  norm  only  if  it  shows  locking  of  at  least  order  f(N) 
in  the  energy  norm. 

(B)  Moreover,  if,  in  addition,  Condition  (a)  holds,  then  5  is  free  from 
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locking  with  respect  to  the  V  norm  iff  it  is  free  with  respect  to  the  energy 
norm.  It  shows  locking  of  order  f(N)  in  the  V  norm  iff  it  shows  locking 
of  order  f(N)  in  the  energy  norm. 

Proof.  Part  (A)  follows  immediately  using  the  fact  that  (2.20)  holds  with  the 
same  Fq(N)  for  both  IMIy  and  11‘llg.  Part  (B)  follows  by  using  Theorems 
2.2  (B)  and  2.3  (B).  □ 


3.  Locking  and  robustness  results  for  the  model  problem. 


In  this  section  we  state  and  prove  some  theoretical  results  for  the  model 
problem,  which  explain  the  computational  observations  presented  in  the 
previous  section. 

We  first  note  that  we  may  write  for  t  €  (0,0]  (0  <  1) 


(3.1) 


bJ(u.v) 


n 


feu 

[dxr 


dv  +  3u  3v_"|  +  fl  _  i  ]  pu 
3x7  +  3xi  dxjjj  |t  J  |3x2 


dx1dx2> 


which  is  of  the  form  (2.22)  with 


a(u,v)  = 


Vu'7v  dx^dx2 

n 


Cu  = 


3u 

3x2' 


Replacing  by  t,  we  have 


t  €  (0,0/1-01  =  (0, b] . 


Clearly  (2.23)  -  (2.24)  are  satisfied  (W  =  L2(Q) ) . 

1  2 

Instead  of  analyzing  the  forms  Bfc  and  Bfc  on  the  same  meshes  as  in 
Section  2,  we  will  now  deal  with  the  analogous  problem  of  considering  only  B^ 
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on  unrotated  and  rotated  meshes  (Figure  3.1).  This  will  simplify  some  nota¬ 
tion.  Although  the  problems  analyzed  now  are  not  totally  Identical  with  the 
ones  addressed  In  Section  2,  they  show  the  same  features  and  are  essentially 
equivalent. 

2 

Let  S  £  R  be  a  triangle  or  paral lelogram  and  let  ^(S)  denote  the 

set  of  all  polynomials  of  total  degree  (degree  In  each  variable)  <  p  if  S 

Is  a  triangle  (paral lelogram) .  We  will  talk  about  spaces  V1^  of  continuous 

piecewise  polynomials  of  degree  5  p  on  the  meshes  in  Figure  3. 1  satisfying 

u|Q  €  7  (S)  for  any  S  in  the  mesh,  for  any  u  e  V**. 
s  p 


Figure  3.1.  Unrotated 

and 

rotated  meshes 

(a) 

Unrotated  triangular  mesh 

(b) 

Unrotated  rectangular  mesh 

(c) 

Rotated  triangular  mesh 

(d) 

Rotated  rectangular  mesh 
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We  will  also  find  it  convenient  to  work  with  a  suitably  scaled  version  of 


(3.1) ,  given  by  the  change  of  variables 

(3.2)  x1  =  Xj,  x2  =  X2V^'  u(x^,  x2)  =  u(x1,x2). 


Under  the  transformation  (3.2)  we  see  that  (3.1)  may  be  written  as 


B, (u, v)  = 


Vt 


Vu*7v  dx^dx2 


a 


with  =  { ( x ^ , x2 ) | I x ^ |  £  1,  I x2 1  £  Vt}  so  that  we  Just  get  Laplace’s  equa¬ 
tion  over  a  very  thin  domain.  The  corresponding  meshes  from  Figure  3. 1  now 
look  as  shown  in  Figure  3.2. 


Figure  3.2.  Scaled  unrotated  and  rotated  meshes. 

(a)  Scaled  unrotated  triangular  mesh  (b)  Scaled  unrotated  rectangular  mesh 

(c)  Scaled  rotated  triangular  mesh  (d)  Scaled  rotated  rectangular  mesh 

The  above  families  of  uniform  meshes  are  all  examples  of  more  general 
quasiuniform  meshes.  We  will  refer  to  triangular  meshes  satisfying  the 
minimum  (respectively,  maximum )  angle  condition  if  the  minimal  angle  of  any 
triangle  in  the  mesh  is  larger  than  c  (respectively,  maximal  angle  is  smal¬ 
ler  than  7i -e),  where  c  is  fixed  for  the  family  of  meshes  (see  [2]). 
Analogously,  we  will  say  that  a  rectangular  mesh  satisfies  the  ratio  condition 
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if  the  ratio  of  the  largest  to  the  smallest  sides  in  each  rectangle  is  bounded 
above  independently  of  the  meshes. 

Let  us  now  look  at  the  scaled  triangular  meshes.  We  observe  that  as 
t— »0,  the  minimum  angle  condition  is  violated  for  both  (a)  and  (c).  How¬ 
ever,  we  observe  a  clear  difference  in  terms  of  the  maximum  angle  condition. 
Mesh  (a)  satisfies  the  maximum  angle  condition  uniformly  with  respect  to  t 
as  h— >0,  while  mesh  (c)  violates  this  condition  as  t— >0.  Similarly,  the 
ratio  condition  for  the  rectangular  mesh  (b)  holds  uniformly  as  t— »0  but 
not  for  the  rotated  mesh  (d).  The  reason  that  meshes  (c)  and  (d)  lead  to 
locking  is  essentially  due  to  the  violation  of  the  maximum  angle  condition  as 
t->0. 


Let  us  define  the  following  weighted  norms  over  12  in  terms  of  the  cor¬ 
responding  norms  over  f2t 


(3.31 


«u«  *  =  ttjIuI  x  . 

H^(n)  tx/*  H*(nt) 


Note  that  ||u|| h1(qj  =  IMIE  t,  the  energy  norm. 

We  will  need  the  following  lemma  about  n-widths. 

Lemma  3.  1.  Let  £2  c  IR1,  1  =  1,2,  be  a  polygonal  domain  (for  1  =  1,  £2  is 

an  Interval ) .  Then 

(3.4)  dN(H1(n),Hk+1(£2)}  *  N_k/i,  k  *1, 

where  the  equivalency  constants  depend  on  k  and  £2  but  are  independent  of 
N. 


For  the  proof,  see  [12]. 

Lemma  3.2.  Let  (v“)  c  H*(£2)  be  subspaces  of  piecewise  polynomials  of  fixed 
degree  p  on  quasi-uniform  triangular  or  rectangular  meshes  which  satisfy  the 
uniform  maximum  angle  condition  (for  triangles)  or  the  ratio  condition  (for 
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rectangles).  Then 


(3.5) 


sup 

ueHk+1(n) 

|u,Hk*1(fl)-1 


inf  || u  -  v |j 
vev  Hx(n) 


>  CN-min(p, k)/2 


where  C  is  independent  of  N,  but  depends  on  p,k,£l  and  the  characteriza¬ 
tion  of  the  mesh  (quasi-uniformity  constant  and  maximum  angle  condition). 


Proof.  For  k  £  p  (3.5)  follows  from  Lemma  3.1.  Let  k  >  p.  In  (3.5)  we 
assumed  that  c  H^fl),  i.e.,  is  the  space  of  conforming  piecewise 

polynomials  of  degree  p.  Let  us  denote  c  L^tfl)  to  be  the  space  of 
piecewise  polynomials  of  degree  p  without  continuity  conditions  between 
elements.  Then  CN  ^  N  £  N  and 


#(N)  =  sup  inf 

u€Hk+1(fl)  veV^ 
l|u||Hk+I(n)11 


[zilu-vll2  I 

Iff  HA  (ST)  J 


1/2 


sup 

ueH1"'!!!) 

|U,Hl“'(n)S1 


inf 

v«v" 


U  -  V 


H!(n) 


where  the  sum  above  is  over  all  single  elements  J  of  the  mesh. 
Hence  it  is  sufficient  to  prove 

4>(N)  2:  CN_p/2. 


-2 

Obviously  also  N  =  N(h)  *  h  «  N.  Consider  now  a  polynomial  u  of  degree 
p+1  (total  (triangles)  or  separate  degrees  (rectangles) ) .  Then  u  e 
H*+*(n).  Hence  if  (3.5)  is  not  true  then 

(3.6)  *(N)  £  f (h)hP 


where  f(h) — >0  as  h — >0.  Since  the  mesh  is  quasi-uniform  with  uniform 
angle  condition  the  following  inverse  inequality  holds  on  every  element  J 
for  any  polynomial  of  degree  at  most  p+1 

||v||  s  Cht_siMI  *.  .  o  S  t  S  s. 

H  (7)  H  (ff) 
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where  C 
have  for 


is  independent  of  y.t.s.h. 


Hence  if  (3.6)  is  satisfied  then  we 


1/2 


if 


u  -  v 


h  HP+ 1 ( 7 ) ' 


<  Cf (h) , 


assuming  that  Hull jjp+1(q)  =  1-  But  since  v^  is  a  polynomial  of  degree  p, 

— >0  as  h— >0,  which  is  a 


ap+1vh 

—  =  0,  i  =  1 , 2  and  hence 


dx 


P+1 

i 


|— r^yll 

3xP  1  L2(«) 


contradiction. 


We  now  consider  the  case  of  the  h-version  on  unrotated  meshes  followed 
by  the  analysis  of  the  rotated  meshes.  For  simplicity  we  will  assume  that 
either  I"..  =  0  or  r_  =  0. 

N  U 

Theorem  3.1.  Consider  the  problem  (2.10)  with  Bt  given  by  (3.1).  Let 
be  the  set  of  piecewise  polynomials  of  degree  ip  (p  £  1)  on  an  unrotated 
mesh  of  the  form  shown  in  Figure  3.1  a, b,  with  mesh  spacing  h  =  h(N).  Let 
Hfc  be  the  space  of  functions  given  by 


H 


t 


{  u  1 1|  u  ||  <  00} 

x(n) 


with  k  1  p.  Let  the  error  measure  Et(u)  =  llull^.  ^  =  ^u^H1(n)‘  Then  t^e 
extension  procedure  5  is  free  from  locking  over  the  region  RQ  =  (0,b]x.y 
and  is  robust  with  uniform  order  g(N)  =  0(N  p/^)  over  this  region. 

Proof.  Let  u  €  H^..  Then  u  €  H^+^(£l^),  with  k  £  p.  Consider  first  the 
case  that  piecewise  linear  functions  are  used  on  a  triangular  mesh.  Then  the 
proof  of  the  sufficiency  of  the  maximum  angle  condition  [2]  easily  gives 


(3.7) 


inf  U  u  -  vf| 


H  (tlt ) 


*  Ch||u|| 


H 2(nt)’ 


where  is  the  image  of  and  (due  to  the  mesh  satisfying  the  maximum 
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angle  condition)  C  is  independent  of  h  and  t. 


Another  interesting  question  is  to  examine  locking  in  the  V  norm  when, 
instead  of  H^+*(£2),  is  taken  to  be  an  unweighted  space.  For  this  case, 

we  expect  once  again  that  there  will  be  no  locking.  We  prove  this  assertion 
for  the  case  p  =  1,  using  the  equivalence  between  finite  difference  and 
finite  element  schemes. 


Theorem  3.2.  Consider  the  problem  (2.10)  with  Bt  given  by  (3.1)  to  be  the 

approximate  problem  corresponding  to  (2.1)  -  (2.2)  with  T..  =  0.  Assume  that 

N 

9  =  (?t>  where  9  =  {/}  is  the  space  of  continuous,  piecewise  linear  func- 

4  - 

tions  on  the  triangular  mesh  shown  in  Figure  3.  1  a.  Let  Ht  =  C  (0) 
(independently  of  t).  Let  the  error  measure  be  E.  (u)  =  ||u||  .  Then  3F 

rtn) 

is  free  of  locking  over  =  (0,b]x./V  and  is  robust  with  uniform  order 
0(N  *^)  there. 


Proof.  Consider  the  mesh  shown  in  Figure  3.  1  a,  with  £1  being  gWen  by 

( —  1 , 1 ) 2  and  the  nodes  being  labeled  A^,  i  =  0,  1 . M,  j  =  0,1,..., M, 

where  M  =  2/h  (Aqq  corresponding  to  the  node  with  coordinates  (-1,-1)). 

N 

Let  Uj.,  the  exact  solution  and  u^.,  the  finite  element  solution  at 

N 

A^j  be  denoted  by  u^j,  u^j,  respectively.  Then  by  the  usual  relation 

between  finite  element  and  finite  difference  solutions  for  a  regular  triangu- 

N 

lar  mesh,  we  have  for  i,j  such  that  A„  €  £1  with  u^  =  0  when  A^  e  T, 
(3.9)  -aJJu£  =  t(-uJ.1>J+2«J  /Vl  j)  *  -  ‘h2fir 

where 


ij 


f*iJ  dxldx2‘ 

£2 


Here,  is  the  standard  "hat  function"  at  A^j.  We  are  given  u  6 

C(4)(£2),  i.e., 


30 


S  K. 


Using  this,  we  may  show 


1  (4)  - 


-  rij +Buh  ■ 


where  =  fCA^)  and  satisfies  |  j  I  S  CK  with  C  independent  of 

u.h  and  t.  Moreover,  for  the  exact  solution  we  have 


(3. 10) 


where  again  |D^j|  £  CK.  Using  (3.9)  -  (3.10),  we  see  that  with  x 
we  have  *  =  0  on  P  and 


(3. 11) 


i  i  4 

-VJ * •  h  Hu 


N 

Ut'V 


with  IH^I  <  CK.  Let 


so  that 


tj(x,  y)  =  1  -  y 


tj.j  =  ( jh) (2-Jh)  >  0. 


Then  =  0  for  j  »  0,M,  and  for  A^  in  the  interior  of  0, 

=  2h2. 

Let  X  =  Ch^Ki?  -  x ■  Then 


(3. 12) 


-aJJX  =  2Ch4K-h4H  2:  Ch4K  >  0. 
t  1  J 


Also,  X  £  0  on  the  boundary  T  of  S3.  Suppose  X.  .  <  0  in  0  for  some 
*  J  ^  J 

A.  «  I\  Then  there  exists  a  minimum  value  of  X  <  0  in  S3.  At  such  a 

A  J  A  J 


value,  for  any  t. 


-a[jx  <  0, 
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contradicting  (3.12).  Hence,  X  >  0  on  0,  i.e., 

A  vl 

*1J  s  “Ay 

2 

Similarly,  taking  £  =  Ch  Ktj  +  x>  we  see  that 

*ij  *  -“Ay 

Noting  that  |t)^|  5  1,  we  have 

IZjjl  5  CKh2. 

This  shows  that  the  difference  between  the  interpolant  of  u^.  and  the  finite 

2 

element  solution  at  the  nodal  points  is  0(h  ),  from  which 

(3.13)  || u  -u^H  ,  <  Ch  =  0(N~1/2). 

1  Hl(0) 

Using  Lemma  3.2.  we  see  that  we  have  the  same  rate  as  (3.5)  with  p  =  1. 

Since  (3.13)  holds  for  all  t,  the  theorem  follows.  □ 

Remark.  In  Theorem  3.2  we  have  used  some  assumptions  only  for  the  sake  of 

simplicity.  For  example,  we  assumed  r  =  Tp,  which  can  be  relaxed.  Further, 

(4)  - 

we  assumed  H^.  =  C  (0)  and  triangular  meshes,  which  can  also  be  relaxed. 

One  could  of  course  consider  the  case  when  H.  =  Hill)  and  E. (u)  = 

V  t 

l|u|lE  £.  To  see  why  this  combination  is  not  appropriate  for  analysis,  consider 
functions  on  Q  that  are  functions  of  the  variable  x^  alone,  i.e.,  4>  e 

Hk(I).  For  such  a  function,  by  Lemma  3.1, 

dM(H1(I),Hk(I))  *  N"(k_1). 

N 

Now  considering  these  functions  as  functions  of  both  x^  and  x^  we  have 
*  €  H^fl)  so  that  we  get  from  the  above 
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Hence,  since 


dM(H1(n),Hk(n),E.  )  >  —  CN~(k_1). 

N  1  Vi 

dN(H1(fi),  H^n),  E^)  does  not  converge  to  zero  uniformly  in  t 
as  N — »oo,  condition  (2.18)  is  violated.  Therefore,  locking  will  exist 
(albeit  in  a  trivial  way)  for  any  choice  of  approximation. 

Let  us  now  consider  the  effect  of  using  rotated  meshes.  Figures  2. 1  - 
2.3  indicate  that  there  is  locking  both  in  the  V  norm  and  the  energy  norm. 
To  explain  this,  let  us  first  prove 


Lemma  3.3.  Let  the  extension  process  5F  be  based  on  the  h-version  using  a 
continuous  piecewise  polynomial  of  fixed  degrees  5  p  on  a  rotated  mesh 
(Figure  3.1  c,d).  Let  u  be  a  function  defined  on  £1  such  that  Cu  =  5u/3x 
=  0  and  u  is  not  a  polynomial  on  Q.  Then  there  exists  a  constant  C  >  0, 
independent  of  N,  such  that 


(3.8) 


inf  || u  -  w|| 

wev 

Cw=0 


inf  || u  -  w||  £  C  for  all  N. 

wef  (£2)  V 

P 

Cw=0 


N 

Proof.  Consider  a  triangular  mesh  of  the  form  shown  in  Figure  3.  1  c.  Let  V 
be  the  space  of  piecewise  polynomials  of  total  degree  £  p  on  this  mesh  and 
let  w  €  satisfy  Cw  =  0.  Then  over  any  element  S7j,  w  must  be  of  the 

form 


w 


i=0 


For  any  two  triangles  7^,3^  with  a  common  side  we  see  that  the  vari¬ 
able  x^  is  never  fixed  along  Hence  the  continuity  of  w  implies  that 

aiJ  =  aik‘ 


Repeating  this  argument,  we  see  that  w  e  .  Since  u  4  JMQ),  we  obtain 
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(3.8).  The  same  argument  works  for  rectangular  meshes.  □ 

We  now  obtain 

Theorem  3.3.  Consider  the  problems  (2.10),  (3.1).  Let  be  the  set  of 

piecewise  polynomials  of  fixed  degree  £  p  (p  £  1)  on  a  rotated  mesh  of  the 
form  shown  in  Figure  3. 1  c  or  d,  with  mesh  spacing  h  =  h(N).  Let  the  error 

measure  be  ^  =  (E*>,  E*(u)  =  || ull r1  ( n)  or  g2  =  ^^t*’  EVU^  =  t’  211101 

let  R={(t,N),  0  <  t  £  b..>  c  R-.  Further,  assume  that  the  quadruple 

N  u 

(3f, ?,S, R)  is  F^-admissible.  Then  the  extension  procedure  ?  shows  complete 
locking  and  is  not  robust  with  any  uniform  order  over  R. 


Proof .  We  can  directly  apply  Theorem  2.2(A)  or  Theorem  2.3(A).  Lemma  3.3 
shows  that  the  condition  (2.32)  is  not  satisfied.  Hence  ?  cannot  be  free 
from  locking  over  R.  The  lack  of  robustness  follows  as  well  from  Lemma 
3.3.  □ 

Suppose  we  now  choose  Ht  to  be  a  space  of  the  form  H^Q),  as  in 
Theorem  3.2.  Then,  since  the  space  of  admissible  functions  in  Theorem  3.3  is 
more  restrictive,  we  again  see  that  there  will  be  locking  in  the  H*(£l)  norm. 
Let  us  mention  here  that  an  alternate  proof  of  Theorem  3.3  may  be  derived  by 
considering  the  necessity  of  the  maximum  angle  condition  (see  [2]). 

Let  us  now  examine  what  happens  when  we  use  the  p-version  instead  of  the 
h-version.  We  have 

Theorem  3.4.  Consider  the  problems  (2.10),  (3.1),  with  =  0.  Let  be 

the  set  of  piecewise  polynomials  of  degree  <  p  =  p(N)  on  any  of  the  meshes 
shown  in  Figure  3.1  (i.e. ,  rotated  or  unrotated),  with  fixed  mesh  size  h 
independent  of  N.  Let  be  the  space  of  functions  given  by 
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Then  the  p-version 


with  k  >  0.  Let  the  error  measure  be  E^(u)  *  ||u||g  ^ 

extension  procedure  5  is  free  from  locking  over  the  region  =  (O.blxA' 

-k/2 

and  is  robust  with  uniform  order  f(N)  =  0(N  )  over  this  region. 

Proof.  Let  u  €  H^+1(Q).  Denote  n  =  <  l x ^  |  Ix^l  <  1>,  ft  =  {Ix^  < 

~  5 

Ix^l  <  £> .  By  the  extension  approach  of  Nikolsky  and  Babich  (see  (11]  Theorem 

~  5 

3.9)  we  extend  u  to  u  on  fl:  for  ”4  <  xi  < 

k+1 

u(iri2)  =  Xxiu['1"(ii+1)i^I,i2] 

1=  1 


with  the  exact  form  of  the  coefficients  being  described  in  [11].  u  is 

„  5  .  .  . 

analogously  extended  on  1  <  x  <  ^  and  then  on  Q  to  get  u(x^,x2),  which 

preserves  the  norm,  i.e., 


>+1(ft) 


5  C ||  u| 


H^+1(n) 


Now  we  use  the  same  approach  as  in  [3]  Lemma  3. 1.  We  transform  ft  onto 

Q  =  {  1^1 1  <  2’  ^2'  <  2^ 

with  the  transformation 


x!  =  4  sin  ?1 


x2  =  ?  sin  C2 


and  let  v(£^,<;2)  =  u(x^,x2).  Then  obviously 


VV+1(Q)  '  CI'GV+1(n)  *C,UV+W 


Now  function  v  can  be  extended  periodically  on  R  =  (|?^|  <  n,  |^2I  <  n> 
and  it  is  symmetric  with  respect  to  the  lines  ^  ~  ^2  =  Wr*t*nS 
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tn 


v  =  Re 


00 


i(JCl+^2) 


and  the  statement  follows  from  (3.10)  and  (3.9).  □ 

Remark.  The  statement  and  proof  of  Theorem  3.4  above  carry  over  identically 
if  instead  of  (3.3),  we  use  the  following  weaker  definition  of  the  norm 
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<<n>’ 


(3. 11) 


u 


I 

0<at^k 


,3aiu„2 


dx 


at"  0 


+  t 


-1 


I 

0Sai+a2^k 

a2£l 


,2. 

a*««ay-° 


In  fact.  Theorem  3.1  also  holds  with  the  above  weaker  definition  -  instead  of 
using  the  maximum  angle  condition,  one  can  establish  it  directly  by  consider¬ 
ing  interpolation  on  the  region  n.  Hence  the  results  in  this  section  are 
valid  if  (3.3)  is  replaced  by  (3.11). 

For  the  case  of  Just  one  square,  the  above  theorem  guarantees  that  the 
rate  of  convergence  is  the  same,  independent  of  t,  for  both  the  rotated  and 
unrotated  case.  This  is  clearly  observed  in  Figure  2.4.  Note  that  our 
function  there  was  u^.  =  sin  e  2  which  satisfies 

uj  €  hJ(B)  V  k. 


Since  this  function  is  very  well  behaved,  we  can  actually  assert  a  stronger 
robustness  result  for  it  (see  Theorem  3.5  below). 

For  the  case  of  more  than  one  square,  Theorem  3.4  asserts  that  the  rate 
of  convergence  is  independent  of  t  both  with  rotated  and  unrotated  meshes. 
However,  in  actual  practice,  while  the  unrotated  case  will  be  insensitive  to 
t,  the  rotated  case  will  not.  This  is  because  it  is  known  from  [4]  that  the 
error  for  the  finite  element  method  using  piecewise  polynomials  of  degree  <  p 
on  a  quasi-uniform  mesh  with  mesh  size  h  is  of  the  form 


(3. 12) 


i  Nn  Ch‘ 

|  u  -  u  ||  £  — 

H  (Q) 


min(k, p) 


u 


Hk+1(n) 


When  the  p-version  is  used,  the  factor  (for  p  large) 

appears  as  a  constant.  In  terms  of  locking,  when  t  is  bounded  away  from  0, 
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one  observes  the  full  effect  of  this  constant.  However,  for  t  close  to  0, 
Lemma  3.3  shows  that  this  constant  is  much  larger.  This  explains  the  data  in 
Table  2.3.  For  a  rigorous  demonstration  of  this  effect,  we  begin  by  proving 
some  lemmas  on  the  approximabl 1 ity  of  smooth  solutions  like  (2.5). 


Lemma  3.4.  Let  f(x)  €  C^d),  I  =  (-1,1),  be  such  that  for  any  k  £  0, 

(3. 13)  Max|f(k)(x) |  £  Adk, 

with  0  <  d  £  1  and  A  Independent  of  k  and  d.  Then 


(3. 14) 


inf  ||f  -v||  <  CAdP+1p(p), 

ve?  ( I ) 

P 


where 


p(p)  = 


2(p+1)(p+l)! 


Moreover,  suppose  that  there  Is  a  B  >  0  such  that  f^p+1^(x)  does  not 


change  sign  on  I  and 
(3. 15) 


minjf ^P+1 ^ (x) |  £  BdP+1. 


Then 

(3.16)  inf  ||f  -  v||  >  CBdP+1p(p), 

ve?  ( I )  U 

P 


where  C  is  independent  of  p.  (Here  ll’llg  -  II  II  ^ •  ) 
Proof.  We  write 


(3. 17) 


r(x)  = 

k=0 


where  L^(x)  is  the  k-th  Legendre  polynomial.  Then  we  have 


.  ,  ,  (-l)k  dk  ...  2.k. 

Lt.  ( x )  =  -n - rr  {(1-x  )  > 

2  k!  dx 
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and  coefficients 


in  (3. 17)  satisfy 


\  = 


2k+l 

2 


1 

f  (x)Ljc(x)dx. 

-1 


Hence  we  have 


(3.18)  lajJ 


2k+l 

2k+1k! 


fl 

I  (l-x2)kf(k)(x)dx|  £  Adk 
-1 


2k+l 

2k+1k! 


1 

( l-x2)k  dx. 

-1 


Denoting 


,1 


H(k)  = 


fl  2  k  (k! )  „k+l 

( 1  X  }  dx  '  (2kVTTV  2 


-1 


we  get 


C1k~1/2  <  K(k)  C2k  1/2 


so  that 


I  a,  |  S 


Ck1/ZAdk 

2klc! 


Now 


‘  k=p+l  p+1 


,2k  1 


(2kk! )2 


CA2  Vf<*l2k,  CA2  fd'|2(p+1) 

((pn).)2^^  '((P+1)!)2W 

since  ^  ^  and  the  sum  is  a  geometric  series.  This  gives  (3.14).  On  the 

other  hand, 

_2  2  2 
ET  £  max  a.  =r— -r. 
p  . .  ,  x  2k+l 

K  k£p+l 

Hence  using  (3.15)  and  (3.18),  we  see  that  (3.16)  holds.  c 


For  any  0  <  d  5  1,  let  us  denote  by  the  set  of  all  functions 
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f(x)  €  C  (I)  such  that  for  each  p  2  1  there  exist  at  most  M  intervals 
l[P\  1  =  (M  independent  of  p),  with  JP  =  I  -  U  I ^  having 

J=1  J 

measure  >  r  >  0  (t  independent  of  p)  and 


mi 
xeJ* 


n  |f(p)(x)|  £  Bdp 


and 

max  | f ^ (x) |  3  AdP 
xel 

with  A,B  independent  of  p.  Then  we  have 

Lemma  3.5.  Let  be  a  sequence  of  quasi-uniform  meshes  on  I  =(-1,1) 

N(h)  h  h 

and  let  I  =  U  I,,  I,  =  (x.,x,  ),  h  =  x ,  -  x ,,  C  h  <  h  S  C  h.  Let 

J  J  J  J+1  J  J+1  J  1  j  2 

{Vp>  denote  a  sequence  of  subspaces  of  H1 ( I )  of  piecewise  polynomials  of 

degree  £  p  on  the  mesh  3rh.  Let  f  e  £ ..  Then 

d 


(3. 19) 


C.dphpp(p-1 )  <  inf  I  f  —  v |  <  C_dHh*p ( p- 1 ) , 

*  _ lrP  it  *  t  r  \  & 


PkP, 


veV: 


H  ( I ) 


where  depend  on  t,A,B,M,  but  are  independent  of  d,h,p  and  where 

the  first  inequality  holds  provided  h  is  sufficiently  small. 

Proof.  Let  g(x)  =  f'(x)  and  let 


(3.20) 


e  =  inf  .  ||g  -  w||  . 

p  1,J  w  €?  .  ( I1?)  L(A 

p-1  J  2  J 


Transforming  1^  to  the  standard  interval  I  we  get  by  applying  Lemma  3.4 


(3.21) 


If  Ij  €  JP,  then  also 


,  S  Cdp*1hp*(1/2)p(p-U. 

r  ^  »  J 
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(3.22) 


'  .  ,  2  CdP*1hP*(1/2)p(p-l). 
*  *  »  J 


Let  now  w(x)  be  the  polynomial  of  degree  p-1  which  attains  the 
inflmum  (3.20).  Then  on  let 


v(x)  =  f(x^)  + 


w(x)dx. 


Because  of  the  construction  of  v(x),  by  expanding  f'(x)  into  Legendre 
polynomials  we  have 

v(xJ+1)  =  f(xJ+1). 

Hence  in  this  manner  we  construct  v(x)  on  the  entire  I  with  v(x)  e  H1(I). 
From  (3.21)  we  have  by  squaring  and  summing  over  j 


|f'(x)  -  v'(x)||L  (I)  <  Cdp+1hpp(p-l). 


For  h  sufficiently  small,  we  have  ^  intervals  contained  in  Jp  with  C 
independent  of  p.  Hence  we  also  get  using  (3.22) 


(x)  -  f '  (x) |tQ  >  CdP+1hPp(p-l). 


Let  us  note  that  the  function  sin  x^  €  and  e  z  e  Hence  the 

function  u*  given  by  (2.5)  satisfies 

(3.23)  «  «t  *  *,  ® 

Note  that  for  sin  x^,  M  =  1  and  (3.19)  will  hold  for  h  £  1.  Also,  if 
h  >  1,  then  (3.19)  holds  for  p  even. 

We  now  prove  the  following  result  for  the  case  when  the  solution  set  ft  = 
(Ht>  where  Hfc  Is  given  by  (3.23).  Here  H^  is  in  fact  a  countably  normed 
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space  and  is  characterized  by  t,A,B,d,T,M. 

Theorem  3.5.  Consider  the  problem  in  Theorem  3.4  with  =  1 1  •  Then 

the  h-p  extension  procedure  ?  using  piecewise  polynomials  of  degree  <  p^ 
on  unrotated  meshes  of  mesh  size  h^,  is  free  from  locking  over  the  region 
Rq  =  (0,b]xjV  and  is  robust  with  uniform  order  f(N)  =  0(h^Np(p^“l ) )  there. 

Proof.  We  consider  for  simplicity  only  the  case  «  a.  Let  ufc  e  so 
that 

ut(xrx2}  *  C(x1)z(x2) 

with  £  e  ty  Then  for  any  y  =  vCx^lwCxg)  e  V**, 

(3.24)  lut-y||g>t  =  IIC'z-v'wIIq  +  IuCz'  -vw'IIq 

^  II?' Z  ~  v'zIIq  +  ||v'x  ~  v,w||q  +  -  Cw'  ||q  +  fllCw'  -vw'IIq 

-  IIS'  -  v'W^Wxwl*  lla:  -  w||2||v'  Iig  + 1<  IICIIqIIx'  -w'll* 

+  lie- v||q||w'  ||q>. 

Suppose  v.w  e  Vp  are  the  best  approximations  to  £,*,  respectively,  in  the 
sense  of  Lemma  3.5,  l.e., 

lie'  -v'ii0  s  chpp( p-i ) 

Hz'  -w'  II 0  S  Ctp/2hpp(p-l). 

Then  (3.24)  gives  for  p  >  1, 

(3.25)  ||ut -y||E>t  ^  ChPp(p-l) 

uniformly  in  t.  This  proves  the  robustness  result.  To  show  that  there  is  no 
locking,  we  must  verify  that  (3.25)  is  the  best  rate  possible.  Now  it  may  be 
easily  shown  that 
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iify  |ut  -  y«!  t  =  lnf „p,ft'vw|lE ,t 

yev  v, weVi 

fl 

*  inf  IIC'Z-v'xll l  =  inf  ||?'  -v'||2||z||2 

veVf.  veV^ 

h  h 

Using  Lenuna  3.5  then  shows  that  (3.25)  is  the  best  rate  possible.  Then  by 
Definition  2. 1,  there  is  no  locking  over  R^.  a 

Remark.  Note  that  the  above  theorem  also  shows  that  there  is  no  locking  for 
the  p-version. 

We  now  look  at  the  case  when  rotated  meshes  are  used.  We  will  distin¬ 
guish  between  two  cases.  For  the  case  that  no  refinement  is  present,  i.e., 
only  one  element  is  used,  we  essentially  get  the  unrotated  case  with  one 

element.  [We  could  understand  this  case  also  by  considering  the  rotated  solu- 
2 

tion  ut  as  in  Section  2  (as  mentioned  in  the  beginning  of  this  section,  this 
is  equivalent).]  We  will  observe  no  locking  in  this  case.  When  refinement  is 
present.  Lemma  3.3  tells  us  that  for  the  limit  problem,  the  mesh  behaves 
essentially  like  one  with  a  single  element.  This  will  show  locking  for  both 
the  p-version  and  h-p  version. 

Theorem  3.6.  Consider  the  problems  (2.10),  (3.1)  (T^  =  a).  Let  be  the 

set  of  piecewise  polynomials  of  degree  p  =  PN  on  a  family  of  rotated  meshes 
of  the  form  shown  in  Figure  3. 1  d,  with  mesh  spacing  h  =  h^.  Let  = 

i?l®i?v^r  and  E^tu)  =  ||u||E  t  or  HuHh1(q)'  Let  ^  consist  of  the 
p-version  (h^  fixed)  or  the  h-p  version.  For  the  p-version  on  a  single 
element,  ?  is  free  from  locking  and  uniformly  robust  with  order  0(p(p^-l)). 
For  a  refined  mesh  with  mesh  size  h^  S  1,  the  extension  process  9  shows 
locking  of  at  least  order  (2/hN)PN  but  is  uniformly  robust  with  order 
p(pN-n. 
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Proof.  For  a  single  element,  we  have  h  =  2.  We  may  either  consider  the  pair 


11  2  2 

u^,  or  ut>  (*,•)-  In  either  case,  Lemma  3.5  may  be  used,  like  in 

the  proof  of  Theorem  3.5,  to  show  that  there  is  no  locking  and  one  gets  uni¬ 
form  robustness  of  order  f ( N )  =  0(p(p^-l)). 

When  refinement  is  present  with  h  S  1,  we  see  that  for  the  limit 
problem,  by  Lemma  3.3,  the  mesh  behaves  like  a  single  element,  i.e.,  one  with 
mesh  size  2.  However,  for  t  bounded  away  from  0,  one  gets  the  same 
convergence  rate  as  in  Theorem  3.5.  This  implies  that 


F0(N) 


inf  ||  u  —  w  ||  v  =  0(hNPNp(pN~l)), 
weV^ 


but  for  Cu  =  0, 


inf  ||u-w||  (F  (N)) 

W€V 

Cw=0 


0(  (r~ )Pn)  . 


Hence  by  (2.34)  we  get  locking  of  at  least  order  (2/hN)PM.  The  method  2 
is  still  uniformly  robust  with  the  rate  using  one  square,  i.e.,  0(p(pN-l)).D 

In  the  above  proof,  we  used  the  necessity  of  the  condition  (2.32)  to 
prove  that  there  was  locking.  We  mention  here  that  for  this  problem.  Condi¬ 
tion  (a)  does  not  hold  in  general  so  that  (2.32)  would  not  be  sufficient 
except  in  special  cases.  For  example,  taking  Hfc  =  span{uj>  gives 


and 


=  sin  Xj 


u, 


0  E,  t 


0(1) 


so  that  if  |MIe  t  £  C||  *  II  t .  then 

'VVi  '  0(1)' 
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violating  (a).  One  could,  of  course,  define  Jf  such  that  we  only  look  at  a 
class  of  problems  for  which  (a)  holds. 

Let  us  summarize  the  various  locking  phenomena  we  have  dealt  with  in  this 
section  in  a  table.  (U  a  Unrotated,  R  a  Rotated. ) 


Extension  Process  E^.  Locking  Robustness 

(uniform  order) 


h-version,  degree  p,  U 

h£+1(A),  kip 

h1 

No 

0(N_p/2) 

( degree  1 ) 

C(4)(«) 

H1 

No 

0(N_p/2) 

h-version,  degree  p,  R 

Hk+1(n),  k  >  0 

h> 

Yes 

None 

Hk+1(fi) ,  k  >  0 

H1 

Yes 

None 

p-version,  U 

H^+1(n),  k  >  0 

hJ.  h1 

No 

-k/P 

0(N  K/^) 

p- vers ion,  R 

H^+1(n),  k  >  0 

<■  h1 

No 

0(N~k/2) 

p  or  h-p  version,  U 

No 

0(hNNp(pN_1) 

p-version,  R,  1  square 

V*Vt 

< 

No 

0(p(pN-l)) 

p  or  h-p  version,  R 

XlxXVt 

< 

Yes 

0(p(pN-l)) 

refined 

We  see  from  the  above  table  the  effect  of  changing  various  choices  of  the 
extension  procedure,  and  E^. 

Let  us  remark  that  our  results  will  be  valid  for  any  rotated  mesh.  We 
have  considered  the  case  that  the  rotation  is  45°  for  convenience,  since 
locking  will  be  observed  most  readily  in  this  case.  However,  Lemma  3.3  is 
true  for  other  angles  of  rotation  as  well. 
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Our  definition  for  locking  in  Section  2  allows  considerable  flexibility 


in  terms  of  the  various  components  etc.,  which  may  depend  on  the 

parameter  t.  We  now  exploit  this  flexibility  to  analyze  an  example  where 
locking  is  caused  by  the  failure  of  an  inverse  estimate  to  hold  when  the 
meshes  are  non  quasi -uni form. 

Let  I  =  [0,1]  and  t  e  (0,^-3.  For  each  N,  let  h  =  N  \  x  =  jh, 

j  =  0,...,N.  Define  5  =  h(l-2t),  c  =  ht  and  let  be  the  sequence  of 

l  r  N  h 

meshes  shown  in  Figure  4.1  consisting  of  nodal  points  x^.x^.x^.  Let  V  =  V 
c  Hq(I)  be  the  space  of  piecewise  linear  functions  on  V  This  defines  the 
family  Given  u  e  Hq ( I ) ,  for  t  >  0,  we  define  u^  ■  u^  e  - 

v£  by  (2. 10)  with 


Figure  4.1.  The  mesh 

i.e.,  u^  is  the  L^  projection  of  u  on  V^.  This  defines  our  extension 
procedure  ?.  We  choose  =  H^(I)  n  Hq(I)  and  E^(u)  =  ||u||^  and  ask 
whether  locking  occurs  as  t— >0.  Note  that  the  only  component  depending  on 
t  is  f  (due  to  the  mesh).  We  let  RQ  =  {(t,N),  0  <  t  £  . 

Theorem  4.1.  Let  {P^.},  H  and  8  be  as  above.  Then  ?  shows  locking  of 
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0(N)  for  any  region  R  =  {(t,N),  0  <  t  £  b.,>  c  R_  and  is  not  robust  with 

N  u 

any  uniform  order. 


Proof.  First  it  is  easy  to  see  that  the  quadruple  (H,  R)  is  F  - 
admissible  with  Fq(N)  *  CN  1  due  to  standard  results  about  the  H1 ( I ) 
projection  on  V^.  For  any  tj  S  t)q  >  0,  the  mesh  is  quasi-uniform  on  R^ 
so  that  the  space  satisfies  the  inverse  assumption.  This  leads  to 

Assumption  A  being  satisfied  so  that  Theorem  2.1  holds  with  ( N )  as  above. 

Let  us  keep  h  fixed  and  let  t— >0.  Then  it  may  be  easily  verified 


h  .  .  h  _  yh 


that  in  the  L 2  norm,  u“— ■■ »u“  e  V^‘  where  vjj  is  a  set  of  discontinuous 


linear  functions  on  E ( J— 1 ) h, Jh ] .  Also  u“  satisfies  the  limit  problem 


(4.  1) 


Jh 

(J-l)h 


h  . 

uQvdx  = 


Jh 

(J-l)h 


uvdx,  Uq(0)  =  Uq(1)  =  0, 


with  v  being  a  linear  function  on  f(j-l)h,jh],  J  =  1 . N. 

Suppose  ?  is  robust  with  some  uniform  order.  Then  for  t  >  0,  h 
fixed, 

llvi^ll  a  S  ||u||  x  -»  || u  -  u£||  x  S  C, 

with  C  independent  of  t.  Hence  there  exists  {u!1  >,  a  convergent 


subsequence  such  that 


h  — h  ,  ..1  f  t  \ 

Utk^U0  ln  H  (I)' 


Since  u^— >Uq  ln  L^,  we  must  have 


u 


But  this  is  a  contradiction,  since  u^  e 
continuous  ln  general,  i.e.,  u^  <t.  Hq(I). 


h 

O’ 

V^,  satisfying  (4.1)  will  not  be 
Hence  ?  is  not  robust  with  any 
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uniform  order.  Also,  because  Fg(N)  *  N  \  5  shows  locking  of  order  N.  □ 

Remark  4.  1.  Note  that  if  we  choose  E(u)  =  ||u||q,  then  3  will  be  free  from 
locking. 


5.  The  clamped  loaded  beam. 

In  this  section  we  look  at  the  robustness  of  some  discretization  schemes 
for  a  one-dimensional  problem  involving  the  Timoshenko  beam.  This  problem 
(along  with  various  extension  processes  for  it)  has  been  analyzed  in  the 
context  of  locking  in  [1],  [9].  Here,  we  present  a  simplified  alternate 
analysis,  using  our  theory,  to  prove  some  of  the  results  in  these  references. 
We  do  this  by  showing  that  Condition  (a)  is  satisfied,  so  that  Theorems  2.2 
(A)  and  (B)  hold.  We  may  therefore  restrict  attention  to  the  limiting  problem 
which  is  essentially  a  one-dimensional  biharmonic  problem. 

The  problem  we  look  at  is 

(5.1)  '*t  +  F(0t“Wt)  =  f  0n  1  =  (0>1} 

(5.2)  F{^t'Wt)/  "  8  on  1 

(5.3)  0t (0)  =  *t(l)  =  wt(0)  *  wt(l)  =  0. 

Here,  wt  and  represent  the  vertical  displacement  and  rotation,  respect¬ 

ively,  of  the  vertical  fibers  of  the  beam,  which  is  subject  to  a  vertical  body 
force  -tg(x)  (f  can  be  related  in  applications  to  dislocations  or  moments). 
The  thickness  Vt  is  assumed  to  lie  in  the  interval  (0,1]. 

We  may  cast  (5.1)  -  (5.3)  into  the  following  variational  form.  Let  V  = 
Hq( I )xHq( I ) .  Then  we  wish  to  find  u^  =  ^t,Wt^  €  V  sat^s^y^n8 
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(5.4) 


Bt(ut,v)  =  b(ut, v)  +  i(Cut,Cv)  =  F(v) 

for  all  v  =  (tfr.z)  e  V,  with  b(u^ ,  v)  =  (^' ,  \p'  ) ,  Cufc  =  0^-w'.  Here  (•,•) 
represents  the  ( I )  inner  product  and  F  is  given  by 

F(v)  =  <f ,0>  +  <g, z>, 

where  <♦,•>  is  the  duality  between  Cl)  and  H  *(I). 

It  is  easy  to  see  [1]  that 

(5.5)  5»utl?  *  llutll|>t  S  (i -|)l“t»v 

with  the  energy  norm  being  defined  as  usual.  The  bilinear  form  can  be  put  in 
the  form  (2.22)  if,  as  in  (3.1),  we  take  instead  of  t,  giving  for 

u  =  (^,w) 

a(u,  u)  *  b(u,u)  +  (Cu.Cu)  2:  ^||u||y 

so  that  (2.23)  -  (2.24)  are  satisfied. 

In  [1]  it  has  been  shown  that  for  f,g  €  H  *(I)  and  0  <  t  S  1,  the 
following  a  priori  estimate  holds  for  k  =  0,1,...,  for  the  problem  (5.4)  (or 
equivalently,  (5.1)  -  (5.3)) 

(5.6)  'VkM*^CutlkSC(lf,k-l+l*lk-l)' 
where  C  is  a  constant  depending  only  on  k  and 


II  ut  H  k+ 1  =  H*tllk+1  + 


t"k+r 


Hence  a  natural  choice  for  the  solution  set  H^  *  H^  is  to  define  Hfc  k+1 

c  H  »  (Hk*1C I )  A  hJ(I))2  by 

(5.?)  Ht  k+1  -  (ull|u|lt  <  .) 

with  the  norm  ||*|lt  =  IMIt  k+1  given  by 
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'u|'t,k+i  •  "u|WEl|Cu|'k- 


Using  (5.1)  -  (5.2),  we  also  have  for  k  =  0,1 . 


(5.8) 


"f"k-!  +  "gllfc-!  SC||ut“t,kT 


From  (5.4),  we  see  that  as  t — >0,  the  Klrchoff  hypothesis 
(5.9)  Cu  =  0 

gets  imposed  so  that  any  uQ  in  the  limiting  set  must  satisfy  (5.9). 
Moreover,  using  (5.6)  with  k  =  0,  we  see  that  if  f,g  are  held  fixed  in 
(5.4)  and  t— »0,  then 


l|Cutl|Q  =  0(t). 

Suppose  now  we  are  given  a  e  t .  We  show  that  Condition  (a) 

holds.  First  it  is  clear  that  |(uilE  t  1  liuHt  k+1  for  k  >  0.  Define  f,g 
by  equations  (5.1)  -  (5.3)  so  that  by  (5.8),  ||f||^_^,  || g|| ^  <  C  and  denote 

by  (u0.C0)  =  ^o,wo’^o^  the  solution  of 


(5. 10) 


-0"  +  =  f  on  t 

?Q  *  8  00  I 

0O  -  W'  =  0  on  I 
0Q(O)  =  0O(1)  =  wQ(0)  =  w0(l). 


(5.  10)  has  the 

(5.  11) 


equivalent  variational  form;  find  (uQ,^0)  €  VxS 
b(uQ,v)  +  (<Q,Cv)  =  f (v)  for  all  v  e  V, 
(CUq.tj)  =  0  for  all  tj  e  S, 


satisfying 


with  S  =  Lgd).  Then  (5.11)  has  a  unique  solution  which  satisfies 
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(5. 12) 


>ub>k+1*«VkSC(>f>k-1  +  >«>lc-l)- 

(5.12)  is  proven  in  [1]  for  k  =  0  and  can  be  obtained  for  higher  k  from 

(5.10)  by  differentiation.  (5.12)  shows  that  uQ  €  HQ  n  Ht  Using  (5.8), 

(5.13)  *  "<o'k  S  C»Vt,k.r 

Now  by  (5.1)  and  (5.10),  we  have 

(5-14)  -<VV"*^‘VV  ■  <0' 

which  corresponds  to  (using  (5.4)  suid  (5.11)) 

(5.15)  b((ut-uQ),v) +^(C(ut-uQ),Cv)  =  (<Q,Cv) 


from  which,  putting  v  =  ut-uQ,  we  get 


(5.16)  ||u  -u, 


,2  .  1 


■  2  -  1 1,  _  „  2  ,  1 . 


l0!lE,t  *  5lVo  +  2F|C(V'V),l,o  s  i|l<ollo  +  fl|ut-uollE,t 


Hence 


'Ut  ”  U0^ E,  t  ~  t^0^0’ 


i.e.  , 
(5. 17) 


lt~'a0*E,t*t  ^O^O’ 


Next,  differentiating  (5.14)  and  using  (5.2),  (5,10)  we  see 


( ■ 0 


so  that 


where  a^ ,  b^ 


(W"  =  V  (W'  =  atX  +  bt 

are  constants  depending  on  t.  Using  (5.17), 


from  which  (using  (5.1)  -  (5.3),  (5.10)  and  (5.14)),  the  following  estimate 
may  be  established 


Using  (5.13),  this  gives 


|ut  ‘  u0*t,k.l  s  1  ll?0*k' 


'“t'Vt.krl  S 


so  that  (a)  is  satisfied.  Hence  we  may  take  to  be  the  set  of  all 

Uq's  constructed  from  ut  e  k  through  the  equations  (5.10). 

We  will  be  interested  in  locking  in  the  ll*l|y  norm.  Then  by  Theorem 
2.2,  a  necessary  and  sufficient  condition  for  the  absence  of  locking  is  that 


for  all  u  a  H^. 


(5. 18) 


inf  || u  -  w||  £  C  Fn(N)  =  infM  ||u-w||v. 

V**  V  0  we/  V 


we 
Cw=0 


We  will  denote  by  (S^>  a  sequence  of  C  subspaces  consisting  of  piecewise 
polynomials  of  degree  <  p  on  a  quasi -uni form  family  of  meshes  on  I. 

Theorem  5.  1.  Consider  the  problem  (5.4).  Let  (V^  =  (S^)^>  =  5^.,  where 
p  >  1  is  kept  fixed  and  the  h  extension  is  used  (h  =  h^).  Let  H  =  {Ht>, 
=  Ht  k+1  be  given  by  (5.7),  with  k  2  p.  Let  the  error  measure  Et(u)  = 
||u||y.  Then  the  h-version  extension  procedure  ?  shows  locking  of  order 
h^1.  It  is  uniformly  robust  with  order  h^  1  for  p  >  1. 

Proof.  We  must  estimate  the  two  inflmums  in  (5.18).  First  we  have  for  u  = 


(5. 19) 


FQ(N)  =  inf  {It-tl.+lw-zl.}  *  0(hj) 
*  Z€Sh 


by  standard  results,  since  <p,  we  H^+^(I),  kip.  Next,  we  estimate  for 


u  =  (0,w)  e  H®  k+1  and  <p  -  w'  =  0, 
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Obviously, 


(5.20)  =  inf  {||w'  -z'  ||  +  ||w-z||  }  *  inf  ||w-z||_. 

^  zssP  1  1  zesP  2 

xl  a 

For  p  *  1,  =  0(1).  For  p  >  2,  standard  spline  theory  [13]  shows  that 

for  kip, 

(5.21)  Aj,  =  0(hP_1). 

(5.19),  (5.21)  together  with  Theorem  2.2  prove  the  result.  □ 

Let  us  now  consider  the  p-version. 

Theorem  5.2.  Consider  the  problem  (5.4).  Let  =  (sj^)2}  =  where  h 

is  fixed  and  the  p-extension  (p  =  pN)  is  used.  Let  fc+1,  k  >  0, 

auid  E^Ju)  =  llu||y.  Then  the  p-version  extension  procedure  ?  is  free  from 
locking  and  uniformly  robust  with  order  p^  . 

Proof.  By  standard  p-version  approximation  theory  we  have,  for  u  =  (0,w)  6 
H® 

o,k+r 

F  (N)  =  inf  <  110  -  ^11  J  +  ||w  -  zll  *  0(p'k) 

+’ZeSh 

since  0,  w  €  H^^d).  Next,  for  u  e  ^+i»  Cu  =  0  so  that  w'  =  <f>  e 
Hk+1(I),  i.e.,  w  €  Hk+2(  I ) .  Hence 

A^  =  inf  |j  w  —  z  II 2  =  0(pjJ(k+2"2))  =  0(p‘k). 

zeSh 

This  proves  the  theorem.  a 

The  results  in  [1],  [9]  also  indicate  what  happens  when  the  measure 
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is  changed.  We  have  summarized  these  results  in  a  table.  The  set  = 

Ht  k+1  in  each  case  (with  k  >  p  for  the  h-version). 

Table  3.1.  Locking  and  robustness  for  different  choices  of  5  and  E^. 


Extension 
Process  SF 

Et 

Locking 

Robustness 
(Uniform  order) 

h-version 

IMIV 

Yes 

None 

P  =  1 

11*11  x 

Yes 

None 

[1] 

llwllj 

Yes 

None 

II  o 

Yes 

None 

I'wIIq 

Yes 

None 

h-version 

llu|lv 

Yes 

OChP'1) 

p  >  2 

11*11! 

Yes 

othP-1) 

111 

llwllj 

No 

(Khj) 

ll^l'o 

Yes 

°(hP) 

"  w"o 

Yes,  experimentally 

O(hP) 

p-version 

llully  • 

No 

°(PNk) 

[9] 

11*11  j 

No 

0(p-k) 

llwllj 

No 

0(p-k) 

■♦■o 

No 

0(p-(k+1)) 

«w|l0 

No 

0(p-<k*u) 

Remark  5.1.  Additional  results  in  (1]  deal  with  the  h-version  of  a  mixed 
formulation  for  which  locking  is  eliminated.  The  reference  [9]  also  analyzes 
an  h-p  extension  procedure  which  is  free  from  locking. 
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6.  Nearly  incompressible  materials. 


The  procedure  outlined  in  the  previous  section  can  be  used  to  analyze 
"Poisson  Locking"  as  well,  which  occurs  in  the  case  of  the  elasticity 
equations  when  the  Poisson  ratio  v  is  close  to  0.5.  More  specif ically, 
consider  the  following  problem 


(6.1) 


-Au  -  grad  dlv  u  =  0  in  0 
v  l-2v  v 


(6.2) 


^  (c,  ,(u.)  +5 

J=1 


U"V  •  -ij  PST  dlv  Vnj  -  «i>  1  s  1  s  2  on  r. 


where  -{ n. j > ^  is  the  unit  outward  normal  to  T  and  the  strain  tensor  (e^) 
is  given  by 


If  a  a  I 

eij(u)  ■  2 [a^j  +  a^i J  • 

It  is  assumed  that 

* 

g*i$ds  =  0 

Jr 

for  any  rigid  motion  R.  Let  denote  t  =  1  2t>  and  u.  =  u  .  As  t— >0, 

v  t  v 

the  constraint 

Cu  =  div  u  =  0 

gets  enforced.  ufc  converges  to  a  limit  u^  which  satisfies 

-Au0  -  grad  =  0  in  (i 

(6.3)  div  uQ  =  0  in  £1 

2 

^c1J(u0)nJ  +  KQni  =  g^  1  £  i  £  2  on  T 

J=1 


which  is  the  Stokes’  equation,  with  being  the  pressure. 
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The  corresponding  variational  forms  for  (6.1)  and  (6.3)  are  respectively 
given  by  (S.4)  and  (5.11)  with  V  =  fij(n),  S  =  <7)  €  L2(Q),  J  7>dx  =  0> , 

2 
IT 
bO 


'n 


«-ii  i 

JJn  i,j-i 


cu(3)Cij(;)dX, 


F(v) 


f  .  . 

g*vds. 


If  the  domain  is  smooth,  we  obtain  all  the  usual  shift  theorems  (as  in 
the  case  of  the  example  in  Section  5).  Also  the  argument  of  equations  (5.15) 

-  (5.16)  carries  over  identically,  so  that  (5.17)  will  again  hold.  In  fact, 
as  we  have  shown  in  [5],  Condition  («)  will  hold  for  this  problem,  even  for 
certain  cases  when  the  domain  is  not  smooth. 

As  a  result,  Theorem  2.2  applies  again  and  we  have  the  following  result 
when  an  appropriate  solution  set  Ht  =  H  is  taken. 

Theorem  6.1.  The  extension  procedure  ?  is  robust  with  respect  to  error 

measure  E^(u)  =  ||u||^,  and  solution  set  H  for  the  variational  form  of  (6.1) 

-  (6.2)  with  uniform  order  g(N)  given  by 

(6.4)  g(N)  =  sup  inf  ||u-w|L. 

Zsf  Zsf 

div  u=0  div  w=0 

It  is  free  from  locking  iff 

(6.5)  g(N)  5  CF  (N)  =  C  sup  inf  ||u-w|| 

SeH5  kv" 

and  shows  locking  of  order  f(N)  iff 

g(N)  *  CFq ( N ) f ( N ) . 

Note  that  in  [6],  (6.5)  was  stated  only  as  a  necessary  condition,  whereas  we 
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show  that  it  is  sufficient  as  well,  which  simplifies  the  subsequent  analysis. 
In  particular,  the  above  result  shows  that  all  the  numerical  schemes  analyzed 
in  [6]  which  were  shown  to  be  robust  for  the  limit  problem  will  also  be  robust 
for  the  sequence  of  problems  given  by  (6.1)  -  (6.2). 

In  15],  we  use  Theorem  6.1  and  some  results  developed  in  Section  2  to 
investigate  the  robustness  of  various  schemes  for  nearly  incompressible 
materials.  In  particular,  we  analyze  the  locking  that  can  result  through  the 
use  of  curved  elements. 
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